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THE MINIMUM INDUCED DRAG OF AEROFOILS. 

13 y Max M. Munk. 


INTRODUCTION. 

The following paper is a dissertation originally presented by the author to the l Diversity 
of Goettingen. It was intended principally for the use of mathematicians and physicists. The 
author is pleased to note that the paper has aroused interest in other circles, to the end that 
the National Advisory Committee for Aeronautics will make it available to a larger circle in 
America. The following introduction has been added in order to first acquaint the reader 
with the essence of the paper. 

In the following development all results are obtained by integrating some simple expressions 
or relations. For our purposes it is sufficient, indeed, to prove the results for a pair of small 
elements. The qualities dealt with are integrable, since, under the assumptions we are allowed 
to make, they can not be affected by integrating. Wo have to consider only the relations 
between any two lifting elements and to add the effects. That is to say, in the process of inte- 
grating each element occurs twice —first, as an element producing an effect, and, second, as an 
element experiencing an effect. In consequence of this the symbols expressing the integration 
look somewhat confusing, and they require so much space in the mathematical expression that 
they are apt to divert the reader’s attention from their real meaning. We have to proceed up 
to* three dimensional problems. TIach clement has to be denoted twice (by a Latin letter and 
by a Greek letter), occurring twice in a different connection. The integral, therefore, is sixfold, 
six symbols of integration standing together and, accordingly, six differentials (always the same) 
standing at the end of the expression, requiring almost the fourth part of the line. The meaning 
of this voluminous group of symbols, however, is not more complicated and not less elementary 
than a single integral or even than a simple addition. 

In section 1 we consider one aerofoil shaped like a straight line and ask how all lifting 
elements, which we assume to be of equal intensity, must be arranged on this line in order to 
offer the least drag. 

If the distribution is, the best one, the drag can not be decreased or increased by transferring 
one lifting element from its old position (a) to some new position (b). For then cither the 
resulting distribution would be improved by this transfer, and therefore was not best before, or 
the transfer of an element from (b) to i(i) would have this effect. Now, the share of one clement 
in the drag is composed of two parts. It takes share in producing a downwash in the neighbor- 
hood of the other lifting elements and, in consequence, a change in their drag. It has itself a 
drag, being situated in the downwash produced by the other elements. 

W (2) 

F,g / 

Considering only two elements, Fig. 1 shows that in the case of the lifting straight line the 
two downwashes, each produced by one element in the neighborhood of the other, are equal. 
For this reason the two drags of the two elements each produced by the other are equal, too, 
and hence the two parts of the entire drag of the wings due to one element. The entire drag 
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produced by one element has twice the value as the drag of that element resulting from the 
downwash in its environs. Hence, the entire drag due to one element is unchanged when the 
element is transferred from one situation to a new one of the same downwash, and the distribu- 
tion is the best only if the downwash is constant over the whole wing. 

In sections 2 to 6 it is shown that the two parts of the drag change by the same value in 
all other cases, too. If the elements are situated in the same transverse plane, the two parts are 
equal. A glance at Fig. 2 shows that the downwash produced by (1) at (2), (3), (4), and (5) 

(2) W 

m 

( 3 ) ( 5 ) 
c=L. 

r,g 2 

is equal. But then it also equals the downwash due to (4), say, produced at (1). This holds 
true even for the component of the downwash in the direction of the lift if the elements are nor- 
mal to each other (Fig. 3.); for this component is proportional x.y/r 3 , according to the symbols 



of the figure. Hence, it is proved for lift of any inclination, horizontal and vertical elements 
being able, by combination, to produce lift in any direction. 

There remains only the question w T hether the two parts of the drag are also equal if the 
elements are situated one behind the other — that is to say, in different longitudinal positions. 
They are not; but their sum is independent of the longitudinal distance apart. To prove 
this, add in Fig. 4 to the lifting element (2) a second inverse lifting element (3) with inverse 



linear longitudinal vortices in the inverse direction. The reader observes that the transverse 
vortices (2) and (3) neutralize each other; the longitudinal linear vortices, however, have the 
same sign, and all four vortices form a pair of vortices running from infinity to infinity. The 
drag, produced by the combination of (1) and this pair, is obviously independent of the longi- 
tudinal positions of (1) and (2). But the added element (3) has not changed the drag, for (1) 
and (3) are situated symmetrically and produce the same mutual downwash. The direction 
of the lift, however, is inverse, and therefore the two drags have the inverse sign, and their sum 
is zero. 

If the two lifting elements are perpendicular to each other (chapter 5), a similar proof can 
be given. 

Sections 6 and 7 contain the conclusions. The condition for a minimum drag does not 
depend upon the longitudinal coordinates, and in order to obtain it the downwash must be 
assumed to be constant at all points in a transverse plane of a corresponding system of aero- 
foils. This is not surprising; the wings act like two dimensional objects accelerating the air 
passing in an infinite transverse plane at a particular moment. Therefore the calculation 
leads to the consideration of the two dimensional flow about the projection of the wings on a 
transverse plane. 
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Section 8 gives the connection between the theory in perfect fluids and the phenomenon 
in true air. It is this connection that allows the application of the results to practical questions. 

1. THE LIFTING STRAIGHT LINE. 

A system of aerofoils moving in an incompressible and frictionless fluid has a dracr (in the 
direction of its motmn) if there is any lift (perpendicular to the direction of its motion). The 
magnitude of this drag depends upon the distribution of the lift over the surface of the aerofoils. 
Although the dimensions of the given system of aerofoils may remain unchanged, the distribu- 
ti°n of the llft can be radically altered by changes in details, such as the aerofoil section or the 
angle of attack. The purpose of the investigation which is given in the following pages is to 
determine (a) the distribution of lift which produces the least drag, and (b) the magnitude of this 
minimum drag. 

Let us first consider a single aerofoil of such dimensions that it may be referred to with 
sufficient exactness as a lifting straight line, which is at right angles to the direction of its flio-ht 
ihe length or span of this line may be denoted by l. Let the line coincide with the horizontal, 
or x axis of a rectangular system of coordinates having its origin at the center of the aerofoil 
lhe density of the lift 


where A, the entire lift from the left end of the wing up to the point x, is generally a function of 
x and may be denoted by f (x ) . Let the velocity of flight be v 0 . 

The modern theory of flight 1 allows the entire drag to be expressed as a definite double 
integral, if certain simplifying assumptions are made. In order to find this integral, it is neces- 
sary to determine the intensity of the longitudinal vortices which run from any lifting element 
to infinity in a direction opposite to the direction of flight. These vortices are generally 
distributed continuously along the whole aerofoil, and their intensity per unit length of the 
aerofoil is 


r ,_ 1 AA' 

v o • p dx ® 

where p is the density of the fluid. Now, for each lifting element dx, we shall calculate the down- 
wash w } which, in accordance with the law of Biot-Savart, is produced at it by all the longi- 
tudinal vortices. A single vortex, beginning at the point x, produces at the point x = f the 
downwash 


dw = 


4 ? rpV 0 


dA f 


Therefore the entire downwash at the point £ is 


1 

i-x 


(3) 


+ T 

1 CdA f i 

“ 4,rpJ dx ■ (4) 

l 

The integration is to be performed along the aerofoil; and the principal value of the integral is 
to be taken at the point x = $. This rule also applies to all of the following integrals. Hence it 
follows that the drag according to the equation 


d W 1IV w 
dx = w ~v. 


is 


"'= 4 4 / fw ■ ( lp 


A'd\ 


1 8ce L - Prandtl, Tragfliigeltheorie, I. Mltteilung. Nachrichten dcr Ges. d. Wiss. zci Gdttingen, 1918 . 


(5) 


( 6 ') 
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or, otherwise expressed, 


i , / 

r 7 + 7 


,r -4i..p J <6 > 

1 i 

~7 “7 

/' here signifies the derivative of/ with respect to x or £. The entire lift is represented by 


A — f f (x) dx 


(7) 


Hence the solution of the problem to determine the best distribution of lift depends upon 
the determination of the function/ so that the double integral 


r 7 +7 


j,=J vw dxd z 




( 8 ) 


l ( 

7 7 


shall have a value as small as possible; while at the same time the value of the simple integral 

+ 2‘ 

J 2 = ft (x)dx = const. (9) 

i 

~ 2‘ 

is fixed. 

The first step towards the solution of this problem is to form the first variation of J 1 


<5 J i — J*j$/ (x) dx [jz I </ + / j ©<« [{/ if 


(10) 


The second integral on the right side of (10) can be reduced to the first. By exchanging the 
symbols x and £ and by partial integration with respect to x , considering/' (£) as the integrable 
factor, there is obtained 




l 

~7 


l 

"7 


/ 

"7 


(11) 


The second member disappears since f—0 at the limits of integration. 2 Further, the right 
hand part, of (11) 

, i 




upon substitution of the new variables x and t = x—t for x and £, is transformed into 


d Cfix-t) 


dx 


J 


dt 


2 If this were not true, there would be infinite velocities at these points. 
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Now 


d 

dx 


2 

/ 


fax-*) //) - f ( 2, 
r d i~ ^ r 

x 2 ^ 9 


h 


f 


(*-*) 


dt 


or, since/ disappears at the limits of integration, 




(x-t) 


dx .1 t 


dt 


which, upon the replacement of the original variables, becomes 

fa* 


so that, finally, 




Substituting this in (10) there finally results 

4 4 * ' 

5./, = 2 J «/(*)<& J /fa? 


•t ' 


( 12 ) 


(13) 


From which the condition for the minimum amount of drag, taking into consideration the 
second condition (9), is 


+, 


JT-** +k "° 


or, when equation (4) is taken into consideration 

w = const. = Wo 


(14) 


(15) 


The necessary condition for the minimum of drag for a lifting straight line is that the down- 
wash produced by the longitudinal vortices be constant along the entire line . 

That this necessary consideration is also sufficient results from the obvious meaning of 
the second variation, which represents the infinitesimal drag produced by the variation of the 
lift if it alone is acting, and therefore it is always greater than zero. 


2. PARALLEL LIFTING ELEMENTS LYING IN A TRANSVERSE PLANE. 

The method just developed may be applied at once to problems of a more general nature. 
If, instead of a single aerofoil, there are several aerofoils in the same straight line perpen- 
dicular to the direction of flight, only the limits of integration arc changed in the development. 
The integration in such cases is to be performed along all of the aerofoils. However, this is 
nonessential for all of the equations and therefore the condition for the minimum drag (equa- 
tion 15) applies to this entire system of aerofoils. 
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Let us now discard the condition that all of the lifting lines are lying in the same straight 
line, but retain, however, the condition that they are parallel to each other, perpendicular to 
the line of flight as before, and that they are all lying in a plane perpendicular to the line of 
flight. Let the height of any lifting line be designated by 2 or Equation (3) transforms into 
a similar one which gives the downwash produced at the point x, 2 by the longitudinal vortex 
beginning on the lifting element at the point f : 


dw= A — dA ' — s > * — 

4 7 rpVo (£ -x) 2 + it-z) 2 

The expression, which must now be a minimum, is 

J t 2) J • /(£, r) (! rz x y- + A -l z y d Z dx 

with the unchanged secondary condition 

J 2 = J/(x, z) dx = const. 

These integrals are to be taken over all of the aerofoils. 

This new problem may be treated in the same manner as the first. 


(3a) 


(8a) 


(9a) 


f “X 

l£-x) 2 + (T-i) 2 


is always to be substituted for It may be shown that this substitution does not 

affect the correctness of equations (10) to (15). Therefore 


w — const. = w 0 (15a) 

is again obtained as the necessary condition for the minimum of the entire drag. 

Finally, this also holds true for the limiting case in which, over a limited portion of the 
transverse plane, the individual aerofoils, like Venetian blinds, lie so closely together that 
they may be considered as a continuous lifting part of a plane. Including all cases which 
have been considered so far, the condition for a minimum of drag can be stated: 

Let the dimensions of a system of aerofoils he given, those in the direction of flight being small 
in comparison with those in other directions . Let the lift he everywhere directed vertically. Under 
these conditions f the downwash produced hy the longitudinal vortices must he uniform at all points 
on the aerofoils in order that there may he a minimum of drag for a given total lift. 


3. THREE DIMENSIONAL PARALLEL LIFTING ELEMENTS. 


The three-dimensional problem may be based upon the two-dimensional one. Let now 
the dimensions in the direction of flight be considerable and let the lifting elements be dis- 
tributed in space in any manner. Let y or v be the coordinates of any point in the direc- 
tion of flight. For the time being, all lifting forces are assumed to be vertical. 

The calculation of the density of drag for this case is somewhat more complicated than in 
the preceding cases. Consideration must be given not only to the longitudinal vortices, which 
are treated as before, but also to the transverse vortices which run perpendicular to the lift at 
any point and to the direction of flight. Their intensity at anv point where there is a lifting 
element is 


F = A ' * - -- 1 =/ ( x , y, z) • - • 

VP w V 0 p 


The density of drag, IF 1 now has two components, W t and W 7 , the first being due to the trans 
verse vortices and the second to the longitudinal vortices. 


100 



MINIMUM INDUCED DRAG OF AEROFOILS. 


For the solution of the present problem only the total drag of all lifting elements 

W=j W'dx 

is to be considered. In the first place it will be shown that the integral of those parts of the 
density resulting from the transverse vortices 

W t = f W/d x 

does not contribute to the total drag. A small element of one transverse vortex of the length dx 
at the point (x, y , z) produces at the point (£, r/ ? f) the downwash 

dw = ^ P v < , 1, ^ Vf(X ’ y ’ Z) - dX (1G) 

where 

r 2 = (J - x) 2 + (t? - y) 2 -f (f - z)\ 

Therefore 

w, = 4 Jp V( ij J / (*, y, 2 ) •/ (i, V, t) v dxd£. (1 7) 

This integration is to be extended over all the aerofoils. It is possible to write this expression 
in such a manner that it holds for a continuous distribution of lift over parts of surfaces or in 
space. This is true, moreover, for most of the expressions in this paper. Now, exchanging the 
variables x, y } z, for £ , y, f, in equation (17) does not change the value of the integral, since the 
symbols for the variables have no influence on the value of a definite integral. On the other 
hand, the factor (rj - y), and therefore the integral also, changes its sign. Hence 

W x = _ W, = 0 (18) 

and, as stated, 

W= W 2 . (19) 

Therefore the entire drag may be calculated without talcing into consideration the transverse 
vortices . 

The method of calculating the effect of the longitudinal vortices can be greatly simplified. 
At the point (f , 77 , f) that part of the density of drag resulting from a longitudinal vortex begin- 
ning at the point (x, y, z) is 

F/-4- f/(£, v, Of {x,y,z)-tdx (20) 

U 0 PJ 

where 

Xf ^ / d r 

f dxf’ reS P- d^ 

and 

* 

i = l r y ~ X ds; f = (£-*)’ + („- s) J + (f - z) 2 . (21) 

V 

The entire drag is 

w= J .w 2 'dx = j f (s , v , r ) f (*, y, 2) Wx. (22) 

Now, in the double integral (22) the variables x, y, z may be exchanged with £, y, f , as before, 
without affecting the value of the definite integral. Partial integration may then be performed 
twice, first with respect to £ and then with respect to x. The substitution results in 

W= 4- f (7 y , 2) f (i, n, 0 fdzdt (23) 

U 0 PJ J 
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^ is obtained from \p upon the exchange of variables. Its value is therefore 

a o 

^ = ds ’ t 2 =(Z-x) 2 +(s-y) 2 + (f-z) J - 

y 

When partially integrating with respect to d £ , the integrable factor is /' (f, 77, f ) 

W 7 = - j* I*/ (I, V, r) ^ i-f (X, y, z) dxd$ 

In the subsequent partial integration with respect to dx } the integrable factor = 

w= ~- v \ p j jf ((> v, f) •/' (*’ v < z ) tdxdt. 

Finally, by addition of (22) and (26), there is obtained 

2 w =^r p j jf (£> 0/ (z. y. 2 ) ('P- 4') dxdt. 


( 24 ) 


(25) 


. d 


dx ^ 


(26) 


(27) 


g = V + tj-s may now be substituted in (24) for the variable of integration s. Then t changes to 
t, and with the exception of the sign the integrand in (21) agrees with the resulting one in (22) 






Subtracting (28) from (21) there results finally 


OO 


(28) 


(29) 


Hence, i/'-f and therefore the entire right side of equation (22) is seen to be independent of the 
longitudinal coordinates y of the lifting elements. 

Therefore the entire resistance of a three-dimensional system of aerofoils with parallel lifting 
elements does not depend upon the longitudinal positions of the lifting elements . 

4. LIFTING ELEMENTS ARRANGED IN ANY DIRECTIONS IN A TRANSVERSE PLANE. 


The problem considered in section 2 can also be generalized in another way. For the present 
the condition that all lifting elements be in one transverse plane may remain. However, they 
need no longer be parallel, and the lift may be due to not only a great number of infinitesimal 
lifts dA but also to similar transverse forces dB. In the first place let the direction of all lifting 
elements be arbitrary, but such that there is a minimum drag, and let this direction be an 
unknown quantity to be determined. 

In the present problem it is desirable to consider a continuous distribution of lift over given 
areas instead of lines. The last case can be deduced from the first at any time by passing to 
the limit. 

Let A* =f(x, z) be the density of the vertical lift per unit area, and B r = F(x, z) the density 
of the lateral force per unit area. The lateral force is considered positive when acting in the 
positive direction of the X-axis. Then the density of the transverse vortices has the com- 


P° nent % oP ■^' and v oV 

density of the transversal 


B\ The density of the longitudinal vortex is the divergence of the 
vortex, or f ^ )• The longitudinal vortices beginning 
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at the point ( 2 , z) therefore produce at the point (£,f) the downwash and the transverse velocity 

(3b) 


1 (dA[_dir 

4?r v 0 p \ dx dz 


j 1 / ibJTL U,l3'\ 7 7 

dw = ^ ( -i — ) dxdz 

2 / 


€— as 

r 2 


7 1 AL4' , 7 

du = i^\dx dz) dxdz 


According to the above, the density of the drag is 


d W= A ,U dxdz + # drds 

t’o Fq 


(3c) 


(5b) 


Witli these symbols there results for the total drag the expression 


w - w p (/// J/' / (U) * ,/ <*wfed#r - f f f f F ' (x ’ 3) F ( ^ )Z A dxdz<m + (3°) 

JJJJ/' (■«,*) f (U) Z ~J dxdziW- (*»*)/«,« 

All of these integrals are to be taken over all of the lifting surfaces. Now the first two 
integrals have forms corresponding to the integral in (8), and therefore there is a possibility of 
substituting (12) for these. A similar relation also holds for the last two integrals. For exam- 
ple, the variation of the third integral is 

5 i // J- f (x,z) F (U) z ~r> dxdzd ^= 

JJJJ [«/' (**)• F ({,f) 8 -=f+f (x,z) SF (f,0 Z ^y x dzd^ (31) 

Now in the first term on the right-hand side the variables x and z may be exchanged with £ 
and It may then be partially integrated with respect to d£, the integrable factor being rf/^(£,f). 
This gives 

f f f f s f ( x > 2 ) ■ F (i»f ) Z A dxdzdids = - J J J J 5 / (f,f) • 4 /.’ ( 2 , 2 ) Lj dxdzd£d$ (32 


This may be partially integrated with respect to dz, the integrable factor being 

d $—z_ d £—2 
df r 2 dz r 2 

JJJJ (X,z) F (t,f) -/ dxdzdm = ~ JjJJV«-»' F ' ix,z) dxdzdm (33) 

Hence the first term of the variation of the third integral of (30) can be transformed into the 
second term of the variation of the fourth integral of this equation. In a similar manner the 
two other terms may be transformed into each other. It is therefore demonstrated that the 
variation of the entire drag may be written 

5 W — 2 j jdf ♦ w * dxdz f 2 ffsF u- dxdz ( 1 3 b ) 
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Two problems of variation can now be stated. In the first place limited parts of the surfaces 
may be at our disposal, over which the vertical lift A and the horizontal transversal force B 
may have any distribution. Only the total lift 

A = JJ/ ( x,z ) dxdz = const. (9b) 

will be given in this case. 

Then 

w = const. =w 0 ) u~o (15b) 

is the condition for the least drag. 

If, however, the lifting parts are similar to lines, there is generally one other condition to 
fulfill. It is then required that the lift disappear everywhere along the direction of the aerofoils. 
That is to say, 

/ sin p — F cos p — o (34) 


where p is the angle of inclination of the aerofoil to the horizontal A"-axis. In order to add the 
new requirement (34) a second Lagrange constant n is introduced. The condition for the least 
drag is now 


and after the elimination of a 


a a 

w;+X t a = o, u— . 

cos P sin p 


w cos P -f u sin p — w Q cos p 


(34a) 


(15c) 


the constant 2 X being replaced by — w 0 , as before. In words: 

If all lifting elements are in one transverse plane , the component of the velocity perpendicular 
to the wings, produced by the longitudinal vortices, must be proportional, at all lifting elements, to 
the cosine of the angle of lateral inclination. 


5. LIFT DISTRIBUTED AND DIRECTED IN ANY MANNER. 

/ 

The results obtained previously can be generalized not only for lifting elements distributed 
in a transverse plane but also for lifting elements distributed in any manner in space. That 
part of the total drag resulting from the transverse vortices is, in the general case 

'■-wCJJJW/ (x, y, z)f (£ V, f) U f dxdydzdidvdt 

(17a) 

' JjJJJJ F (x, y, z) F (£, t), f) dxdydzd^drjdfj 

Both terms have the same form as the integral in (17). The demonstration for (17) therefore 
applies to both. In the general case also the total drag can be calculated from the longitudinal 
vortices without taking into consideration the transverse vortices. 


In this as in (20) , 


W= JJj J J J J jfo, y, z) f (£, v, f) dxdydzdtdydt 
+ J J J [JJ F-(x,y, z) F (£, 0 dxdydzdidr\dS 

- / J J J J V ’ ^ F dxdydzd£dvd{ 

- / J / J / J” F 0,y,z)f (£, v, f) ^2 dxdydzdtd-od£ J 


00 

‘■‘-U 


t-X 


ds : 


-£JV* 


t- = (t-x ) 2 + (r,-sy + ({-zy 


(22a) 
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The first two terms in (22a) have the same form as the right-liand side of (22), and the same 
conclusions are therefore valid for each. It can he proved directly for (22a) as for (22) that 

each of the two double integrals is independent of the longitudinal coordinates of the lifting 

elements. This proof can now be extended over the last two integrals of equation (22a). 

The third integral, after changing the variables, becomes 

/// J7 // ( x > V’ x ) F ' (f> V, f) <Pi dxdydzd£dr)d<; = (35) 

/ ////// (£> r !’ f) F' (x, V, z) dbhdtdxdydz 

where 


= “ 4 ^. X ds; t 2 = (£ - x) 1 + (s - y) 2 + (f - 2 ) 2 

V 

Now, let /^' be chosen as the integrable factor and be partially integrated with respect to z. 

JJJJJ// (x, y, 2 ) I ’ ($, V, f) iidxdydzdzdndf^ (3(3) 

JJJJJ// & * « * (*. 2 ) £ iAdvh'dxdydz. 


As in the previous cases, the second integral to be expected vanishes since / as well as F 
disappear at the limits of the integration. Next ^ ^ is chosen as the integrable factor 

and partially integrated with respect to x . By \j/ 2} by analogy, is meant 



y 


JJJJJ j/ ^ x ’ y< F) F' ($, V, f) 'b,dxdydzdtdT)d$ = (37) 

JJ J JJJ/' (S> V, f) ^ (z, ?/, 2 ) iJtdrjdpixdydz. 


Now \p 2 may be transformed, the variable s in the defining equation being replaced by v + $• 

The result is that 


^2 = 4“ J" - 1 3 ” ds. < 2 = (£ - x) 2 + (7; - s) 2 + (f - 2) 2 . 

1/ 

It is seen that the integrand agrees with that of the defining integral \p 2 . Therefore, and since 
the right-hand side of (37) 'contains the same function under the double integral as the fourth 
turn in (22a), this fourth term can be combined with the transformed third member. This 
gives 

J JJ JJJ-/ ( ' x ’ V> F) F’ (£, j;, f) ^ idxdydzd^d-qd^ + (38) 

J/JJ/J*’ <* * 2 )/' (£, V, f ) • $ 2 dxdydzd%dridt; = 

JJJJJJ^ 2/> 2 )/' (£. f) Hi - W dxdydzd^d-qdl 

where 

+ OO 

— 00 

Fi — 'l'z an( l therefore the two sides of (38) are independent of y. This is therefore demon- 
strated for the whole right-hand side of (22a). 

201G7— 23 25 


105 



REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 


In general it can therefore be said: 

The total resistance is always independent oj the longitudinal coordinates oj the lifting elements. 

And further: 

The most favorable distribution oj the lift, with reference to the total drag, occurs when this is 
also the case for the projection oj the lifting elements on a transverse plane . 

That is to say, all of the lifting elements are projected on a plane perpendicular to the 
direction of flight, and any element so obtained has a lift equal to the sum of the lifts of all lifting 
elements projected onto it. 

6. DETERMINATION OF THE SOLUTIONS. 

The previous demonstrations show that the investigation for the distribution of lift which 
causes the least drag is reduced to the solution of the problem for sj^stems of aerofoils which are 
situated in a plane perpendicular to the direction of flight. In addition, the condition for least 
drag (15c), which becomes the condition of uniform downwash (15) if the lift is vertical, leads 
to a problem which has often been investigated in the theory of two-dimensional flow with a 
logarithmic potential. The flow produced within the lifting transverse plane by the longitudinal 
vortices originating in it is, indeed, of this type. Each such vortex produces a distribution of 
velocity such as is produced by a two-dimensional vortex of half its intensity, and the whole 
distribution of velocity is obtained by adding the distributions produced by the longitudinal 
vortices. The potential flow sought is determined by the condition of (15c). Let it be com- 
bined with the flow of constant vertical upward motion w=-a> 0 . The resulting flow satisfies 
the condition at the boundaries 

w cos 6 + (x sin j3 0 (39) 

and there results, for the case of lifting lines: 

The two dimension at potential flow is of the type that encircles the lifting lines , and at a gnat 
distance the velocity is directed upwards and has the value w — — w 0 . s 

Within lifting* surfaces the velocity is zero according to the condition (15b), and the fluid 
therefore flows around the contour. 

The intensity of the longitudinal vortices at any point is twice the rotation of the two 
dimensional flow. In the case of the lifting lines, therefore, the density of the longitudinal 
vortices is double the discontinuity of velocity from one side to the other. The intensity of the 
transversal vortices is determined by integrating the longitudinal vortices along the aerofoils 
and therefore equals jLwice the difference of the velocity-integral produced on the two sides of 
the aerofoil. Now the integral of the velocity produced is identical with the potential and 
hence it appears: 

The density of the lift perpendicular to the lifting line is proportional to the discontinuity of 
potential — and has the value 

JA'* + B'> = 2v 0 p(<pi-‘P 1 ) (40) 

Hence the total lift obtained by integrating over all aerofoils is 

A = 2v 0 pf (<p 2 - Vl )dx (41) 

Sometimes a transformation of this equation is useful. In order to obtain it, suppose that 
all of the lifting lines are divided into small parts. Then, on the two ends of each lifting element 
there begin tv. o in /"rse longitudinal vortices, the effect of which on a distant point is that of a 
double vortex. Their velocity-potential <p and their stream function \p may be combined in the 
comjdex function \ p+i<p, and, not considering the existence of a parallel flow, which is without 
any importance in the calculation, this complex function has the form for a lifting line, 


d{yp-\rif) = 


dA + idB 

z-z 0 
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where 2 represents x -f- iy and 2 0 ---=x 0 + iy 0 , x o and y 0 being the coordinates of tlie lifting elements of 
the line. For a lift distributed over areas a similar equation can be formed. The integration of 


(42) gives 


\p -f i<p = 


CdA + idB 

J Z ~ Zq 


(42a) 


Now the residuum of the integrand at infinity is dA -f idB and therefore the residuum of the 
integral is A-\~iR. Therefore the expression can be written. 


A = 2i\pR [Res (\p + iip) ] 


(41a) 


where the last part means the real part of the residuum of \p + i<p at infinity. In the most im- 
portant case of horizontal aerofoils the residuum itself is real and can be used directly to calcu- 
late the lift. The density of drag at any point is proportional to the perpendical component of 

the density of lift and is — A', from which results 1F= U °. A. Making use of (41) one obtains 

v 0 % 


W- 


AC 





(43) 


m__ A 2 ^0 

2v^p R [Res (^ + ?V)] 


(43a) 


The integral in the denominator of (43) represents an area characteristic of the system of 
aerofoils investigated. Frequently the easiest method of calculation is to assume from the 
beginning the velocity ?r 0 at infinity to be unity. 

The case of the lift continuously distributed over single parts of areas is derived from the 
preceding one by passing to the limit. Since the vertical velocity w disappears at all points in 
the lifting surfaces, the velocity is zero at all points and the rotation vanishes. 

Therefore , in the case of the most favorable distribution of lift, all of the longitudinal vortices 
from the continuously lifting areas begin at the boundaries of the areas . 

Equations (43) and (43a) remain. The distribution of lift is indeterminate to a certain 
extent. On the other hand, it is possible to connect the points of the contour having the same 
potential <p by strips of any form, and it is only necessary that the lift be always perpendicular 
to the strip and its density have a constant value along the whole strip. According to equation 
(40) this equals the difference of the potential at the contour between the two borders of the 
strip. Worthy of note is the special case in which all of the strips run along the contour, thus 
coining again to the case of lifting lines. It. appears that: 

Closed lines have the same mi minium of drag as the enclosed areas when continuously loaded. 
Especially 7 important are those symmetrical contours which are cut by horizontal lines in 
only two points. With such the limitation to vertical lift does not involve an increase of the 
minimum drag. For this case it appears that: 

The density of the vertical lift 'per unit area must be proportional to the vertical component of 
the velocity of the two-dimensional flow at the point of the contour of the same height z . It is 


dA 

dF 


— 2v 0 p 


d<p 

dz 


The corresponding density of drag is 


dW_ 

dF 



(44) 


(45) 
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7. EXAMPLES OF CALCULATIONS. 

Examples of calculation of the previous demonstrations can he based on any calculated 
two-dimensional potential How around parts of lines or areas. The simplest flow of the first 
kind is that around a single horizontal line. It leads to the problem investigated at the begin- 
ning of this paper. 

In this case the potential is the real part of where p denotes x + iz. Ihe lifting 

line joins the two points 2=0, x=-l and 2 = 0, x= +/, and has the length 2. The velocity at 
infinity is w = t. The discontinuity of potential along the lifting 'line is ~ V, = 2\ 1 —x*. Ihe 
density of lift is distributed according to the same law, therefore if plotted over the span the 
density of lift would be represented by the half of an ellipse. 


The minimum drag is 


W=A 2 


1 1 

Vo 2 p/ 2 4tt 


(4G) 


if, instead of the value 2, the span had the general value b , the minimum drag would be 


W- 


1 1 
v 2 p/ 2 7 rid 


(47) 


This same result has been obtained by Prof. Prandtl by another method. 3 
The simplest example for a lifting vertical area is the circle. Let its center coincide with 
the origin of the svstem of coordinates. Then the potential of the flow around this circle is 


<p = 


(48) 


where r=M f x*+z 2 . At infinity W 0 = l. Under the condition of and according to equation (40) 
the density of lift is 

<49) 

This results in a constant density of lift of A' = 2. Therefore the drag is 


A 2 


1 

2 i; 0 2 p 


j j 2d xdz 


, 1 
Vp/ 2 - Stt 


(50) 


The double integral is to be taken over the circle. If the general case for the diameter 
equal to L) be considered, then the least drag is 


Vf=A* 


l 

Up/-.' 


1 

l) 2 2 tt 


(51) 


Hence in respect to the minimum drag the circle is equivalent to a lifting line having 
a length y2 times the diameter. 

A lifting circular line would have the same minimum drag as the circular area. 

This result was also obtained by Prof. Prandtl by another method. 4 A reduction of the 
original problem of variation to the two-dimensional flow sometimes enables a survey of the 
result to be made without calculation. For instance, lot a third aerofoil be added between the 
two aerofoils of a biplane having a small gap. (The gap may be about one-sixth of the 
span.) Then, in order to find the most favorable distribution of lift, the double line about 
which the How occurs is to be replaced by three lifting lines. Now, in the region of the middle 
lifting line the velocity is small, even before this line is introduced. Therefore the discontinuity 
of the potential along the middle line is very much smaller than that along the others. Hence 
it results that the middle aerofoil of a triplane should lift less than the other two. 

* Fi F st communication concerning t nin Zeitschrift fur Flugtechnik uml Motorl. iyi4. S. '239, in a note by Betz. 

< See Technischc Berichle der Fiugzeugmeisterei Bd. II Heft 3. 
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8. PROCEDURE FOR THE CASE OF FLUIDS WITH SMALL VISCOSITY. 

The preceding results do not apply so much to the calculation of the most favorable distribu- 
tion of lift as to the calculation of the least drag. For it appears, and the results are checked 
by calculation, that even considerable variations from the condition of most favorable distribu- 
tion of lift do not increase the drag to any great extent. Usually the minimum drag can be 
considered as the real drag of the system of aerofoils and in order to allow for the effect of 
friction of the air it is sufficient to make an addition. This addition depends chiefly upon the 
aerofoil section; it also depends, omitting the Reynolds Number, only upon the area of the 
wings and on the dynamical pressure. It is independent of the dimensions of the system of 
wings themselves. It may be useful to have a name for that part of the density of drag, inde- 
pendent of the friction of the air, which results from the theory developed in this paper. It is 
called the “induced drag.” Generally it is not the drag itself but an absolute coefficient which 
is considered. This coefficient is defined by 

m 

where It* is the drag previously denoted by W f q is the dynamical pressure vJ.p/2, and F is 
the total area of the wings. Equation (43) can now be written 


7 r(k- b) 2 


where c a is the lift coefficient ^ f corresponding to c w . The greatest horizontal span b of the 

system of wings perpendicular to the direction of flight is arbitrarily chosen as a length char- 
acteristic of the proportions of the system, k is a factor characteristic of the system of aero- 
foils and has, according to the preceding, the value. 


1 | V 2 "' { r 

’ b 2 J U ' 0 


It has a special physical significance. 

Under the same conditions a single aerofoil with a span of k tunes the maximum span of a 
system of aerofoils has the same induced min imum resistance as the system. 

9. REFINEMENT OF THE THEORY. 

The demonstrations given rest on the assumption that the velocities produced by the 
vortices are small in comparison with the velocity of flight. The next assumption, more ac- 
curate, would he that only powers higher than the first power could be neglected. 

In this case the solutions just found for lifting elements in a transverse plane can be con- 
sidered as the first step towards the calculation of more exact solutions. The following steps 

must be taken: The exact density of drag is = ” where v is the horizontal velocity 

produced at the lifting elements by the transverse vortices. It can he calculated exactly 
enough from the first approximation. Now, the condition of least drag is 

w- cos (3 + n sin j3 = ?r 0 cos 0 ^ 1 + - ^ ( 1 5d) 

and the flow of potential, according to this condition at the boundary, is to be found. Compared 
with the first approximation the density below is in general somewhat increased and the den- 
sity above is somewhat decreased. The minimum drag changes only by quantities of the 
second order. 
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If the lifting elements are distributed in three dimensions a similar refinement can easily 
be found. In this case there is to be taken into consideration a second factor which always 
comes in if the differences of the longitudinal coordinates of the lifting elements arc consider- 
able. The direction of the longitudinal vortices do not agree exactly with the direction of 
flight, but they coincide with the direction of the velocity of the fluid around the aerofoil. 
They are therefore somewhat inclined downwards. A better approximation is obtained by 
projecting the lifting elements not in the direction of flight but in a direction slightly inclined 

'2.ir 

upwards from the rear to the front. This inclination is about ^ _ n - . Except for this, the method 
of calculation remains unchanged. 
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